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a = limit of L\Ja(L) at r = 2 calculated from Eq. (12). The FP action data 
scales within its uncertainty and appears to agree with the zero-lattice spacing 
extrapolation. 

We have described a systematic program for constructing fixed point actions 
for SU(3) gauge theory, and illustrated it using a particular scale two blocking 
transformation. An approximation to the FP action shows scaling within our 
small statistical errors beginning at aT c = 1/2, as compared to aT c < 1/8 for 
the Wilson action, at a cost of a factor of about 5-7 in computational speed per 
site update. The outstanding problem now is to extend the program to include 
fermions. The free fermion case, the Schwinger model and the large N limit of 
the Gross-Neveu model have been considered in this context [15]. 
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Figure 3: Scaling test for the Wilson action (crosses) and FP action (squares); T c 
is defined in infinite spatial volumes. The diamonds are from zero-temperature 
simulations by Ref. 14 and the diamond and burst at a = are extrapolations. 



dent of the bare coupling (resolution) in the scaling (continuum) limit. Any 
variation of G is due to lattice artifacts. Figure 3 shows G(L) with r = 2.0 as 
a function of aT c for the Wilson and for the 8 parameter FP action. The inner 
and outer error bars on the FP points show their statistical uncertainty, and 
the combined uncertainty from statistics and in f3 c . The Wilson action shows 
a scaling violation of about 10% between N t = 2 and 6. No scaling violation 
above the statistical errors is seen for the FP action. 

No published simulations with the Wilson action at small lattice spacing 
precisely match our scaling test using the torelon mass on fixed physical volumes. 
The ones which come closest are those of Boyd et. al. [14]. They present 
measurements of the quantity y/a '/T c at the critical couplings for Nf = 4, 6, 8, 
and 12, for which the string tension has been computed via fits to Wilson loops 
on large (zero temperature) lattices. The string tension from Wilson loops is 
an upper bound on <r{L); <r{L) is reduced from a by the zero-point fluctuation 
term which for large L in a string model is 

a(L) = a-^ + .... (12) 

The diamonds in Fig. 3 show 2y/a/T c from Ref. [14]. The authors of that 
work present an extrapolation to a = 0, which we also display along with the 
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procedure has considerable uncertainty in it: we are attempting to construct a 
practical approximation to a FP action subject to the constraint that it be easy 
to use. It is not obvious to us what set of coarse configurations constitutes a 
"reasonable" data set to use in the construction, nor how closely the action must 
approximate a true FP action in order to perform well in simulations. These 
problems are open questions. 

For numerical simulation we need a parametrization which is simple and 
approximates the FP action sufficiently well on those configurations which are 
typical in the actual calculation. In order to find such a parametrization we 
fitted S FP (V) with powers of the traces of only 2 loops. We took the plaquette 
and the twisted 6 link operator (x,y,z,-x,-y,-z). We kept up to four powers of 
(3 — ReTW) n . The 8 parameters are listed in Table 1. 

Table 1: Couplings of the few-parameter FP action for the RG transformation. 



operator 




C'2 


C3 


c 4 


Cplaq 


.523 


.0021 


.0053 


.0167 


Cd — link 


.0597 


.0054 


.0051 


-.0006 



We simulated this action with a mixed Metropolis-overrelaxation algorithm 
acting on SU(2) subgroups, in complete analogy with a standard Wilson action 
code. No special optimization was done. The program runs about a factor of 
5-7 slower than a highly optimized Wilson code. 

Scaling means that all physical dimensional quantities show the same func- 
tional dependence on the gauge coupling. It is different from asymptotic scal- 
ing where in addition to scaling we require that this functional dependence be 
described by the 2-loop A parameter. Asymptotic scaling properties can be 
significantly improved by a suitable non-perturbative redefinition of the bare 
coupling [12, 13]. However, here we are not concerned with asymptotic scaling. 
Scaling properties can only be changed by modifying the lattice action. 

We used the critical temperature T c to set the physical scale for our scaling 
test. We determined the critical coupling constant f3 c (Nt) for N t = 2,3,4 and 
6 and fixed the lattice spacing at these coupling values as a = l/(T c Nt). 

Next we measured physical observables at couplings f3 c (Nt) in fixed, finite 
physical volumes. This way we can avoid infinite volume extrapolations. For 
one scaling test we consider the quantity G = L\J u(L) where u(L) is the string 
tension on an L 3 volume computed from the exponential fall-off of the Polyakov 
line correlator (torelon mass). As we vary the lattice coupling we also vary the 
lattice spatial size so that the physical volume is kept fixed at V = (r/T c ) 3 , 
where r is some conveniently chosen aspect ratio. The quantity G is indepen- 
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Figure 2: The quantity rV(r) constructed for the Wilson action. 

the fixed point action. The minimization of Eq. (3) gives the fine configuration 
{U } for this {V} 

S'(V)= minu(S (U)+T(U,V)) (10) 
S (U ) + T(U ,V). 

If So is close to the fixed point action S FP , than Uo is close to {U FP } and up 
to quadratic corrections 

S FP (V)= S FP (Uo) + T(Uo,V) 

= S FP (Uo) + (S'(V)-So(Uo)). (11) 

This equation can be used to calculate the value of the fixed point action on the 
course configuration {V} if S FP is known on the fine configuration {Uo}. 

We generated about 400 configurations on 2 4 and 4 4 coarse lattices using 
the Wilson action, with coupling [3 ranging from 5.0 to 7.0 and solved Eq. (11), 
parameterizing the FP action as a sum of powers of traces of loops. The fitting 
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2. One can prove that the FP correlator has no power law corrections, only 
exponentially vanishing ~ exp(— r/a) ones. 
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Figure 1: The quantity rV(r) constructed from the correlators of the FP 
Polyakov loops and using the FP action. In the continuum rV(r) = 1/4-7T = 
.0796. 

We have constructed FP actions for two types of RG transformations. The 
one we have studied most extensively is based on a scale factor 2 blocking [11]. 
It has tunable parameters, which we adjust to minimize the range of the action. 
We did this construction in two steps. First, we solved the RG equations for 
smooth gauge configurations. In this limit the action can be taken to be an 
arbitrary quadratic function of the vector potential. The solution of the RG 
equations can be formulated analytically and the computer is only used to do 
Fourier mode sums. Here we do the tuning of parameters for locality. To 
go beyond leading order we solved the fixed point equation Eq. (3) numerically. 
The procedure is the following: We create an arbitrary SU(3) configuration {V}. 
That serves as a coarse configuration. We then minimize the RHS of Eq. (3) 
numerically. As an intermediate step we approximate S FP with a simpler action 
So- We choose So such that it is sufficiently simple for minimization but close to 
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g?i small, c? 2 = c?3 = 
RG step 



= 0, i. e. the classical fixed point action, then after one 



d[ = di + aiVdf + 0(d 3 ) (7) 
4 = a^dl + 0(d 3 ), (8) 

or 

^S FP 1 -WT(S FP + a <$g*R 2 + a <$g*R 3 + . . .) (9) 

9 g 2 + a\{g A 

so that the action itself is changed only in order g 4 to all orders in the latticed 
spacing a. Note that variations on the Symanzik[8, 9] program yield actions 
which are only 0(a 2 ) or 0(g 2 a 2 ) improved. 

An explicit perturbative calculation with the FP action of the non-linear 
cr-model in d = 2 [10] confirms this picture. In ref. [10] the mass gap m(L) was 
calculated in a finite periodic box of size L . The prediction of the standard action 
is contaminated by power like cut-off" effects (a 2 /L 2 )",n = 1,2, .... With the 
FP action these type of cut-off effects are reduced to the level of the numerical 
precision (dominated by the uncertainties in the couplings of the quartic part 
of the action). The coefficient of the leading a 2 j L 2 cut-off effect is reduced by 
4 orders of magnitude relative to that of the standard action. 

The action on the renormalized trajectory is perfect — it has only physical 
excitations. In a Green's function only physical states come in as intermediate 
states independent of the operator used. In spite of that Green's functions (of 
fields, currents, etc.) generally show cut-off effects. Consider, for example a free 
scalar field and imagine constructing the FP action by making many blocking 
steps, beginning on a very fine lattice. The FP lattice action is described by 
a lattice field which is the average of the continuum field over a hypercube. 
Therefore, the two-point function of the lattice field in the FP action is equal 
to the two-point function of the continuum field averaged over the hypercubes 
around the the two lattice points in question. This averaging brings in a trivial 
rotation symmetry breaking. The two-point function of the field in the FP action 
is not rotationally invariant, although only physical states propagate. The same 
result is true for the potential derived from Polyakov loop correlators: the naive 
correlator has power law corrections to the true potential. 

These power law corrections can be eliminated by constructing fixed point 
operators — operators which the naive operators block into under repeated 
RG transformations. We have constructed the FP Polyakov loop for smooth 
gauge configurations (quadratic action) and used its correlator to compute the 
perturbative [3 — ► oo potential V(r). We display the classical potential from 
the FP action using its FP operator, in Fig. 1 and contrast it to the poten- 
tial computed from the Wilson action using the naive operator, shown in Fig. 
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Let us denote the couplings of f3S(U) by /3, 02,03, . . .. The action can be 
represented by a point in the infinite dimensional coupling constant space. A RG 
transformation introduces blocked links on a coarse lattice and integrates out 
the original lattice variables, thus increasing the lattice spacing while leaving 
the physical long distance predictions unchanged. The blocked link variable 
V pins) ESU(N) is coupled to a local average of the original link variables. The 
new action is defined as 



where T(U, V) is the blocking kernel. 

Under repeated real space renormalization group transformations the action 
moves in this multidimensional space. On the critical surface it flows to a fixed 
point. For an asymptotically free theory the critical surface is at /3 — ► 00. In 
this limit Eq. (2) reduces to a saddle point problem giving 



where S FP is the FP action. The FP action is constructed by solving Eq. (3). 
This action is classically perfect: its spectrum exhibits no lattice artifacts. It 
also has scale invariant instanton solutions [6]. 

In addition to being classically perfect, the FP action is also one-loop per- 
fect. Moving away from the fixed point along the RT, the only effect of the 
perturbative g 2 corrections is to make the marginal operator S FP (U) weakly 
relevant, i.e. the coupling /3 = 2N/g 2 begins to move, but the form of the action 
on the renormalized trajectory remains unchanged through 0(g 2 ): 



under a RG transformation, to 1-loop level. Here /?' = /3 — A/3, where A/3 is 
fixed by the first universal coefficient of the /3-function. An equation analogous 
to Eq. (4) was presented in ref. [7] without discussion. 

The formal RG argument is very simple. Parameterize the theory in terms of 
one marginal operator R\, its coupling d\ , and many irrelevant operators R n (U), 
n = 2, . . . and their associated couplings d n . Under an RG transformation 





S FP (V) = mm (S FP (U) + T(U, V)) , 



(3) 



pS FP (U) p'S FP (V), 



(4) 



d' 



'k 



^kdk + a 



fd l d J +0(d 3 ). 



(5) 



where Ai = 1 and A n < 1 for n > 1. If we start with the action 



5" 



,FP 



Ri, 



(6) 
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Most of the nonperturbative studies of quantum field theory used today 
involve introducing a discrete space-time lattice as an ultraviolet regulator. 
Cut-off independent continuum quantities are obtained in the limit of zero lat- 
tice spacing. However the lattice introduces artifacts like rotational symmetry 
breaking. The main problem for numerical lattice calculations is the control of 
finite lattice spacing effects. 

The influence of lattice artifacts on numerical results can depend strongly on 
the specific regularization or lattice action. Most numerical calculations of QCD 
to date use the Wilson plaquette action. This action is the simplest one which 
can be used but simulations done with it do not show scaling unless the lattice 
spacing is less then 0.1 fm and, consequently, nowadays typical pure gauge or 
quenched calculations are done on very large spatial volumes [1]. 

Two of us [2] proposed using a lattice action which is completely free of lattice 
artifacts. That such a "perfect" action exists follows from Wilson's renormal- 
ization group theory [3]. The fixed point (FP) of a renormalization group (RG) 
transformation and the renormalized trajectory (RT) form a perfect quantum 
action. The FP action itself reproduces all the important properties of the 
continuum classical action (it is a classical perfect action). Although at finite 
coupling the FP action is not "quantum perfect" it is expected to be a very 
good first approximation. For the 2-dimensional non-linear a model an action 
parametrized by about 20 parameters was found that showed no finite lattice 
spacing effects even at correlation length £ ~ 3 lattice spacings, at a negligible 
computational overhead compared to the gain obtained by the great reduction 
of lattice artifacts. 

This note is a summary of work we have presented in two papers [4, 5] 
constructing and testing a FP action for lattice SU (3) pure gauge theory. The 
goal of this program is to find lattice actions which can be used in simulations 
and which show improved scaling properties compared to the Wilson action. 
Along the way we discovered some remarkable properties of FP actions, whose 
relevance may extend beyond lattice calculations. 

The partition function of an SU(N) gauge theory defined on a hypercubic 
lattice is 



where f3S{U) is some lattice representation of the continuum action. It is a func- 
tion of the products of link variables U^(n) = e lA ^ n "> GSU(N) along arbitrary 
closed loops. The normalization is fixed such that on smooth configurations 
the action takes the standard continuum form f3S(U) — ► / d A xTi {F^F^), 
where f3 = 2N/g 2 . 
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Abstract 

We summarize our recent work on the construction and properties 
of fixed point (FP) actions for lattice SU (3) pure gauge theory. These 
actions have scale invariant instanton solutions and their spectrum is exact 
through 1-loop, i.e. in their physical predictions there are no a n nor g 2 a n 
cut-off effects for any n. We present a few-parameter approximation to a 
classical FP action which is valid for short correlation lengths. We pe rform 
a scaling test of the action by computing the quantity G = L^J cr{L), 
where the string tension cr(L) is measured from the torelon mass fi = 
La(L), on lattices of fixed physical volume and varying lattice spacing a. 
While the Wilson action shows scaling violations of about ten per cent, 
the approximate fixed point action scales within the statistical errors for 
1/2 > aT c . 
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